Observation of Cs Rydberg atom macrodimers 
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We report the observation of cold Cs Rydberg atom molecules bound at internuclear separations of 
R ~3-9 fira. The bound states result from avoided crossings between Rydberg atom pair interaction 
potentials in an applied electric field. The molecular states can be modified by changing the applied 
electric field. The molecules are observed by mapping the radial separation of the two Rydberg 
atoms as a function of time delay between excitation and detection using the Coulomb repulsion of 
the ions after pulsed field ionization. Measurements were performed for 63D + 65D, 64D + 66D, 
65D + 67D, and 66D + 68D pairs. The experiment is in good agreement with calculations of the 
pair interactions for these states. 

PACS numbers: 34.50.Cx,32.80.Ec,82.20.Bc 



Frozen Rydberg gases 0, Q have been the subject of 
intense research recently. The construction of fast quan- 
tum gates and single photon sources using dipole block- 
ade [1, 0, d, 0, 0JR H, [Toll . the study of cold Rydberg 
atom molecules [ill, Il2l I la, Il4 |. and the investigation 
of many body physics [2J are central motivations for 
this work. Cold Rydberg atom molecules are exciting 
because of the interesting properties that these objects 
possess. Molecules formed by two cold Rydberg atoms 
are called macrodimers since the atoms are bound at dis- 
tances > 1 fim. It has been suggested that due to their 
delicate nature, macrodimers can be used to study vac- 
uum fluctuations, quenching in ultracold collisions [13] 
and Rydberg atom interactions including their controlla- 
bility with applied electric fields. Prior experiments have 
been unable to unambiguously confirm that these unique 
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states of matter exist as bound states 
the first experimental observation of bound macrodimers 
which have one of the largest, if not the largest, molecular 
bond observed to date. 

The macrodimers that we observe result from avoided 
crossings between Rydberg atom pair interaction po- 
tentials in an applied electric field, e These 
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FIG. 1: Integrated atomic ion yield spectra and pair poten- 
tials for 65D + Q7D with e = 190 mV cm - . The excitation 
laser intensity is ~ 500 W cm" 2 . All fine structure and Q are 
plotted. f2 = mji+m-,2 is the projection of the angular mo- 
mentum on R. The dashed line highlights one of the features 
studied in this paper. 



macrodimers are unique because they feature a quasi- 
continuum of bound states, are found at extremely large 
internuclear separation, R ~3-9 /im, and are formed by 
the interplay between Van der Waals interactions and the 
Stark effect resulting from e. e can be used to stabilize, 
destroy or modify the potential well that gives rise to 
the bound states. A spectrum and calculated potentials 
[ItJ for the 65D + 67D pair with e = 190mVcm~ x , are 
shown in Fig. [TJ The potentials have prominent wells 
at R ^3-9 /xm. These potential wells support hundreds 
of bound states [j~3rj with maximum energy spacings of 
~100kHz. The lifetimes of the molecules arc limited by 
the radiative and blackbody decay of the atoms [3] • 



DETECTING MACRODIMERS 

Macrodimers are difficult to detect. Pulsed field ion- 
ization (PFI) breaks the molecule apart and complicates 
time-of-flight (TOF) mass spectrometry. Using selective 
field ionization to ionize only one of the Rydberg atoms 
forming the molecule creates an ion close to the remain- 
ing Rydberg atom. The electric field of the ion can ionize 
the atom. Field ionizing both atoms creates two ions that 
repel each other, causing the complex to dissociate. All 
these possibilities lead to the detection of singly ionized 
Cs atoms at the same average TOF. Detection of the 
molecular ro-vibrational spectrum is difficult because it 
is a quasi-continuum. For the state shown in Fig. [TJ the 
rotational constant is ~ 6 Hz while the maximum clas- 
sical angular momentum at the equilibrium separation, 
Re, is Imax = /ivR e /h w 310. [i is the reduced mass 
and v is the thermal velocity. Given that the radiative 
decay time leads to a natural line width of ~ 10 kHz for 
66-D and taking into account the allowed rotational tran- 
sitions due to the large impact parameter, each vibra- 
tional line will have a width of ~ 580 kHz. This width is 
much greater than the vibrational energy spacings. Con- 
sequently, it is unlikely that the ro-vibrational states of a 
macrodimer can be resolved unless deeper and narrower 
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potential wells can be located. 

We use the recoil velocity distribution to detect the Cs 
macrodimers in our experiment [lij ]. A basic char- 
acteristic of a bound state is that the atoms remain in 
close proximity for a long time compared to the vibra- 
tional period. The signature of molecule formation in 
our experiment is based on this principle. We measure 
the recoil velocity of the ion fragments after PFI as a 
function of the delay, r, between excitation and PFI. If 
two Rydberg atoms are close together and pulsed field 
ionized, the ions created will Coulomb repel each other. 
The amount of Coulomb repulsion, which is reflected in 
the recoil velocity of the ions, is a measure of the R at 
which the pair was ionized. Measuring the recoil velocity 
of the ions as a function of r maps the distribution of ra- 
dial separations of the two Rydberg atoms as a function 
of time. For a bound state, the distribution of R will 
remain relatively constant as a function of r. For a dis- 
sociative process, the distribution of R will change with 
r. The recoil velocity can be measured by recording the 
TOF distribution of the ions that were created by PFI of 
the molecule at a detector in a TOF experiment [l8j |. 

In our experiment, we do not expect to observe oscil- 
lations of the Coulomb repulsion due to molecular ro- 
vibrational wave-packet motion. The vibrational periods 
are < 10 /is. Experimental limitations such as long term 
laser stability and excitation pulse average away any os- 
cillations that may be present. Future experiments de- 
signed to observe wave-packet motion will be interesting 
since macrodimers exist at the interface between quan- 
tum and classical dynamics. 



EXPERIMENTAL METHODS AND RESULTS 

The macrodimers are created in a Cs magneto-optic 
trap (MOT). The MOT contains 10 7 atoms at a de nsity 
of p ~ 2 x 10 10 cm -3 at a temperature of ~ 124 pK [llf . 
The background pressure is ~ 5 x 10~ 10 Torr. The MOT 
is formed between two electric field plates, 26.7 cm above 
a microchannel plate (MCP) detector. The experimental 
apparatus is described in more detail in (l8j . 

The Cs macrodimers are excited by a 2-photon scheme 
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6P 3/2 + 2hv (n - l)D + (n + l)D, (1) 



where n = 64, 65, 66, and 67. The MOT trapping light is 
shut off during the time when the molecules are excited 
and detected. The MOT repumping light is on during 
the experiment. The atoms in the 6P 3 / 2 (F = 5) state 
are excited from the 65 , 1/ / 2 (P = 4) state using a 1mm 
spot size beam derived from the trapping laser. This 
beam is crossed at 22.5° with a beam from a dye laser 
(Coherent 699) operating at A = 508 — 509 nm that drives 
the 2-photon transition in Eq. [TJ The dye laser is focused 
to a ~ 25 ± 10 pm spot size that was measured with a 
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FIG. 2: The rate dependence on excitation laser intensity. 
Solid triangles (a) denote 63D+65D molecular state, open 
triangles (A) 64D+66D, open circles (o) the 65D+67D, and 
closed squares (■) 66D+68D. Intensities and rates have been 
normalized to show the quadratic dependence. The solid line 
is a quadratic fit to the data. 



CCD camera. The dye laser polarization is parallel to 
the TOF axis. The crossed beam geometry helps reduce 
background counts. The intensity of the dye laser and the 
trapping light used to excite the atoms to 6P 3 / 2 (P = 5) 
is 500 Wcm -2 and 0.5 Wcm~ 2 respectively. 

Bound or free Rydberg atoms are detected in the ex- 
periment using PFI. The experiment begins when the 
trapping light is switched off. 2.5 ps after this light is 
switched off, the Rydberg excitation beams are turned 
on for 5 ps. After a variable delay, r, a PFI pulse is 
applied for 2 ps at an amplitude of 76 V cm" 1 to ionize 
the Rydberg atoms that were excited and push the ions 
to the MCP. The rise time of the electric field pulse is 
~ 80ns. A 153mVcm _1 ps^ 1 , bps duration ramp elec- 
tric field is applied immediately preceding the PFI pulse 
to clear any stray ions from the excitation region. The 
atoms move < 500 nm due to thermal motion during this 
time. The ionized atoms have a TOF distribution cen- 
tered at ~ 27 ps. The pulses from the MCP are processed 
by a constant fraction discriminator and accumulated by 
a multichannel analyzer. A constant voltage is applied 
to the upper field plate during excitation to generate e. 
To measure TOF velocity distributions, the sequence was 
repeated at a rate of 1 kHz which resulted in count rates 
of - 10 Hz. 

Under the experimental conditions described, many 
molecular resonances are observed for Rydberg atom 
pairs at different e [T3,[3|. In this experiment, we mea- 
sured molecular resonances for 63D + 65D, 64D + 66D, 
65D + 67D and 66D + 68D pairs. We focused on these 
states because our calculations indicated that bound 
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FIG. 3: Comparison of TOF distributions for the 66D atomic 
case (blue) and the 65D+67D molecular case (red) at r = 
18 /is. The horizontal axis is centered on the atomic TOF. 
These distributions measure the particle velocities or kinetic 
energies [18J . The inset shows the full width at half maximum 
of a Gaussian fit to each of the different distributions at dif- 
ferent time delays. The atomic data was taken by tuning the 
laser to the atomic state under conditions of single particle 
counting. 



states existed for certain e and these states satisfied the 
necessary experimental requirements. The experimental 
considerations are a compromise between an acceptable 
oscillator strength and the elimination of background 
counts from Rydberg atoms. The oscillator strength for 
the molecular excitation is predominantly due to mix- 
ing via multipole interactions and electric fields with 
the nearest atomic state which tends to decrease with 
increasing energy separation. The asymptotic atomic 
states are far detuned from the excitation laser fre- 
quency so without the mixing one expects little transition 
strength to the molecular states. However, the noise in 
the measurement increases with n as the molecular reso- 
nance moves closer to the atomic line wing where some 
atomic Rydberg excitation occurs. 

The excitation rate for each molecular resonance was 
first measured as a function of excitation laser intensity 
and found to be quadratic indicating a 2-photon excita- 
tion process, Fig. [2] These observations, our prior exper- 
iments [l6[ and the fact that the gas is dilute are evidence 
that the spectral features are from processes involving 
two Rydberg atoms. The fact that the spectral features 
are associated with molecular excitations does not nec- 
essarily indicate the presence of macrodimers. Many res- 
onances result from photo-initiated collisions where the 
atoms are not excited to bound states [l6j |. 

Two sample TOF distributions at r = 18 /zs are shown 
in Fig. [3J The inset of Fig. [3] shows the full width at half 
maximum of a Gaussian fit to each of the distributions for 
r = 18 — 38/is. The narrow distribution was obtained by 



tuning the excitation laser to the 66D atomic resonance 
under conditions for single atom counting. The broad 
distribution was acquired by tuning the excitation laser 
to the feature circled in Fig. [TJ The difference between 
the two distributions in Fig. [3] is a result of the dynam- 
ics taking place when single atoms versus pairs of atoms 
are excited. Rydberg atoms expand from the excitation 
volume at the thermal velocity, ~ 8.9 cm s _1 for Cs. Cs 
atom pairs expand at the vector sum of their thermal ve- 
locity, collisional recoil velocity and the velocity acquired 
from Coulomb repulsion after PFI. The molecular signal 
is much broader than the atomic signal, characteristic of 
the Coulomb repulsion between the ions created at close 
range, R < 10 /mi, by PFI. Although the amplitude of 
the traces shown in Fig. [3] is similar, the count rates for 
the atomic signal are 3-4 orders of magnitude larger for 
identical excitation laser intensities. 

The TOF distribution of a molecule will only be broad- 
ened by the thermal velocity and the velocity obtained 
from Coulomb repulsion because the collisional recoil ve- 
locity is zero. If the time between PFI and the average 
arrival time at the detector is not changed, the Coulomb 
repulsion between the atoms forming a molecule will re- 
main constant as a function of t since the distribution 
of R remains constant. The TOF distribution under 
these conditions measures the distance between the ex- 
cited atoms as a function of r. This description is valid as 
long as there is no observable wave-packet motion which 
would lead to an oscillating signal. 

Pairs of atoms undergoing a collision have increas- 
ing internuclear separation, i?, with increasing t so the 
Coulomb repulsion will decrease as r increases for fixed 
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FIG. 4: This figure shows a plot of kinetic energy vs. r for 
a collision between two 89D atoms. The line is the energy 
expected for a pair of Cs atoms starting at R = 2.8 ^m and 
moving uniformly apart at a velocity of 17 cms -1 per atom 
for a time r before ionization as measured for this collision in 
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FIG. 5: Recoil kinetic energy vs. delay for 63D + 65D with e = 224 mV cm -1 , 64D + 66D with e = 205 mV cm" 1 , 65D + 67D 
with e = 190mVcm _1 and 66D + 68D with e = 158mVcm _1 . Spheres denote the data taken for the molecular resonances 
and the dotted black line shows a calculation of the expected results of a binary collision with thermal velocity as described in 
the text. 



TOF. An example of this behavior is shown in Fig. 2] for 
the collision between 89D atoms studied in |l6[. These 
dynamics will result in a Coulomb repulsion broadened 
TOF distribution at short r that becomes smaller as R in- 
creases due to the exit channel velocity of the fragments. 
Asymptotically, R > 15 /im, when two atoms leaving a 
collision are ionized, the Coulomb repulsion is negligi- 
ble and the velocity is determined by the exit channel 
velocity and sample temperature alone [l6| . Collisions 
where the atoms are excited at long range and are al- 
lowed to evolve on attractive molecular potentials can 
also be studied in this way but this is a subject of future 
work. 

To observe macrodimers, TOF distributions were ac- 
quired for t — 18 — 38 [is. Figure O shows the data for 
the 4 states used in the experiments. The times were 
chosen to be short enough to avoid decay of one of the 
molecular partners, but long enough to measure the de- 
pendence on t so that a collision could be distinguished 
from a macrodimer. The lifetime of a ~ 60-D state includ- 
ing blackbody decay is ~ 100 /is 2(| 21 1. Qualitatively, 
the kinetic energy distribution is not changing with r in- 
dicating that the spectral feature that the laser is tuned 
to is a macrodimer. 

A linear fit was applied to the data for each molecular 
pair shown in Fig. [5] to determine the minimum colli- 
sion velocity that would give a similar TOF distribution. 
The velocity allowed by the experimental error in the 
measured distributions is ± 1cm s _1 for 63D + 65-D, 
-2.7 ± 2.1cm s" 1 for 64L> + 66D, 2.2 ± 2.7cm s" 1 for 
65D + 67 D and 3.4 ± 4 cm s" 1 for 66L> + 68D. All errors 
quoted are standard deviations. The collision velocities 
are all zero within the resolution of our spectrometer, 
~ 2.5 cms^ 1 , for all pairs li|, li| . The noise in the mea- 



surement increases with n because the molecular states 
are getting energetically closer to the nearest atomic line. 

The data is also consistent with our calculations of the 
interatomic potential wells. Notice that the Coulomb 
repulsion decreases with increasing n. This is because 
the potential wells shift to larger R as n increases. The 
fact that the 64D + 66D and 65D + 67D macrodimers 
have approximately the same recoil kinetic energy is also 
predicted by our calculations. The R e of the wells closest 
to the spectral lines studied at the e quoted in Fig. [5] for 
63L> + 65L>, 64D + 66L>, 65L> + 67£>, and 66L> + 68D 
are approximately 3.3 /im, 3.5 fmx, 3.6 //m and 3.9 /im, 
respectively. The R e are a measure of the potential well 
shape. It should be pointed out, however, that most of 
the excitation probability is at large R. This is due to the 
Frank-Condon factors and the number of available pairs 
at particular R at our densities. These factors favor near 
dissociation states. The macrodimers spend most of their 
time at R — 4 — 9 /im at the outer turning points of the 
potential wells. We also note that there is not a simple 
n scaling law dependence to the well shape because this 
depends on e as well as n. 

We compared the measured TOF distributions to sim- 
ulations of the ion trajectories in our apparatus. Each 
simulation includes the spectrometer geometry, e, and 
the PFI pulse. R for each macrodimer was chosen uni- 
formly from a range based on the calculated potential 
wells. Molecules or atoms were placed at random in the 
excitation volume according to the excitation laser inten- 
sity distribution. The simulated TOF distribution width 
was in good agreement with the molecular signal. The 
atomic measurement yields a ~ 50 % larger TOF distri- 
bution most likely due to the alignment of the dye laser 
focal spot to the MOT and the uncertainty in determin- 
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ing the laser spot size. The Coulomb repulsion for our 
range of r is at least 3 times larger than the thermal 
broadening so the spot size uncertainty does not prevent 
the measurement. The expected behavior shown in Fig. 
H] as a solid line is also in good agreement with the data 
further verifying the accuracy of the simulations. 

We compared the experimental data on the 
macrodimers to the expected behavior of an ultra- 
low energy, 5.22 MHz, binary collision. The recoil 
velocity of the atom pair was taken to be equal to the 
thermal velocity of the cold gas and the initial R was 
set equal to that determined by the Coulomb repulsion 
measured in the experiment at the minimum r. The 
calculation is shown in Fig. [5] as a dotted line. The data 
from the experiment clearly indicates that the recoil 
velocity of the pair of atoms excited in the experiment is 
much less than the thermal velocity of the gas. Note that 
larger collision velocities cause the recoil kinetic energy 
from the Coulomb repulsion to decrease faster. This 
experimental result and the fact that our calculations 
predict bound states gives us confidence that we have 
observed bound macrodimers [l6| . The measurement for 
4 different states and the fact that the potential energy 
calculations predict potential wells that support bound 
states makes it highly unlikely that the laser, 1.5 MHz 
spectral bandwidth, could selectively excite pairs with 
relative velocities of < 3.4 ms . 



modifications to the apparatus. We plan to study this 
effect in the future. 

To conclude, we have demonstrated the existence of 
bound molecular states of Cs Rydberg atom pairs by 
mapping the R distribution of the atoms using the 
Coulomb repulsion of the ions after PFI. Our experi- 
ments are in good agreement with Monte Carlo calcu- 
lations using no adjustable parameters of the TOF dis- 
tributions and calculations of the pair interactions of the 
Rydberg atoms. These macrodimers are stabilized by e 
and result from avoided crossings, e can be used to con- 
trol the formation and characteristics of the molecules. 
Macrodimers like these will be observable in other species 
of cold Rydberg gases. With the higher densities that can 
be achieved in an optical dipole trap, imaging studies can 
be carried out to investigate the alignment effect noted 
here. Now that macrodimers have been observed, their 
production can be optimized and they can be trapped 
so these sensitive, novel states of matter can be used for 
new experiments. 
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DISCUSSION 

The potential wells probed in our experiment depend 
strongly on e. The resonances vanish at larger e because 
the potentials have become repulsive. At e = OmVcm -1 , 
our calculations do not predict any of the prominent wells 
shown in Fig. [1] and the spectral features are not ob- 
served. In the cases explored in this paper, the molecu- 
lar states are shifted as e is increased until they interact 
to form avoided crossings. As e is increased further, the 
states are shifted through each other so the bound states 
created by the avoided crossings are destroyed. The elec- 
tric fields were chosen for our experiments so that the po- 
tentials wells were broad and deep so that the molecules 
were most easily observed. 

A closer look at the molecular signal in Fig. [^indicates 
that there may be an alignment effect. The symmetric 
maxima about the peak center were observed for all the 
data. If there were no alignment effect, the TOF dis- 
tribution would look flat topped with a smoothing from 
the thermal velocity [lj|. We measured the TOF dis- 
tribution for the molecular resonance with the dye laser 
polarization perpendicular to the TOF axis to explore 
this feature. We observed no difference in the width or 
shape of the distribution. We conclude that if alignment 
of the molecule is present, it is due to e. We currently 
have no way of verifying this conclusion without major 
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